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Abstract 

The critical value of the atom-field coupling strength for a finite number of atoms is deter- 
mined by means of both, semiclassical and exact solutions. In the semiclassical approach we use 
a variational procedure with coherent and symmetry-adapted states, while for the exact quantum 
solution the concept of fidelity is employed. These procedures allow for the determination of the 
phase transitions in the model, and coincide in the thermodynamic limit. For the three cases men- 
tioned above, universal parametric curves are obtained for the expectation values of both the first 
quadrature of the electromagnetic field, and the atomic relative population, as implicit functions 
of the atom-field coupling parameter, valid for the ground- and first-excited states. 

PACS numbers: 42.50.Ct, 03.65.Fd, 64.70.Tg 



I. INTRODUCTION 



A many-body system (e.g. a cold 2-level atomic cloud) interacting with a 1-mode radi- 
ation field inside an optical cavity in the dipolar approximation is described by the Dicke 
model pQ . An important feature of this model is the presence of a phase transition from the 
normal to the superradiant behavior [2j. This is a collective effect involving all A atoms in 
the sample, where the decay rate is proportional to A 2 instead of A, the expected result for 
independent atom emission. While this transition has been much debated in the literature, 
mainly due to the fact that it requires a very strong atom-field coupling, recent experimental 
results indicate that it may actually be observed in situations in which atomic transitions are 
replaced by transitions between momentum states [3, 4J. This has brought about a renewed 
interest in the Dicke model, mainly for its phase transitions but also because it can present 
multi-partite entanglement, and can be realized in systems more widely than in the original 
cavity QED case (cf. e.g. [5] and references therein). 

In this work we obtain the phase transition of the Dicke model for a finite number A of 
atoms, via 3 different methods: i) through a numerical diagonalization of the Hamiltonian 
and the use of the fidelity between neighbouring states; ii) through variational test states that 
are a direct product of coherent Heisenberg-Weyl (ii/"W(l))-states (for the electromagnetic 
field), and S r ?7(2)-states (for the atomic field) P,[7j; and hi) through the use of projection 
operators on the coherent states in (ii) to obtain states which obey the parity symmetry 
in the total excitation number, present in the Hamiltonian. All these procedures coincide 
in the thermodynamic limit A oo. We show that, for a finite number of atoms, the 
symmetry- adapted states (iii) constitute a much better approximation to the exact quantum 
ground- and first excited-states. (Needless to say, the symmetry is also exploited in the 
numerical diagonalization of the Hamiltonian.) 

For the three cases mentioned above, universal parametric curves are obtained for the 
expectation values of both the first quadrature of the electromagnetic field, and the atomic 
relative population, as functions of the atom-field coupling parameter, valid for the ground- 
and first-excited states. 
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II. PHASE TRANSITIONS AND UNIVERSAL BEHAVIOR 



The Dicke Hamiltonian involves the collective interaction of TV two-level atoms, with en- 
ergy separation fudAi with a one-mode radiation field of frequency Cjf ) in the long wavelength 
limit. It has the form 

H D = a ] a + uj a J z + -J= (a f + a) ( J+ + J_) , (1) 

where uja = uja/^f and 7, the coupling parameter between the matter and field, are given 
in units of the frequency of the field, and we have taken h — 1. The operators a\ a denote 
the one-mode creation and annihilation photon operators; J z the atomic relative population 
operator; and j± the atomic transition operators. 

The energy surface is found by taking the expectation value of the Dicke Hamiltonian 
with respect to the tensorial product of coherent states for the HW(1) and SU(2) groups 
\a) <g> |C) , given by EE] 

\a) = exp(-|a| 2 /2) ^~j=M > 

where the parameters a and £ are complex numbers. 
Straightforward calculation leads to [TO] 

= ^ (p 2 + <? 2 ) ~ i u A cos 6> + 2 ^7 g sin 6> cos cj) . (2) 

In this expression we use the harmonic oscillator realization for the field variables and the 
stereographic projection for the angular momentum parameters, 
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a = -^(q + ip) , C = e"^ tan - , (3) 

where (q,p) correspond to the expectation values of the quadratures of the field, and (0, (j)) 
determine a point on the Bloch sphere. 

Using the Ritz variational principle one finds the best variational approximation to the 
ground state energy of the system and its corresponding eigenstate. The minima and degen- 
erate critical points are obtained by means of the catastrophe formalism. For the coherent 
states we find that when 7^ = glu/4 the critical points degenerate; for this critical value of 
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the field-matter coupling, a phase transition takes place from the normal to the superradiant 
behavior of the atoms pj] . The critical points which minimize T~L are given by 



0r = O- 



q c = 0, 



p c = 0, for |7| < 7 C , 



(4) 



6 C = arccos(7 c /7) 2 , q c = y/l- (7 c /7) 4 cos C , p c = , for \j\ > j c , 

where (f) c = 0, 7r. Fig. [l] shows q c /\/JV and # c as a function of 7, for (f) c = 0, ecu = 1 and 
TV = 20 atoms. The phase transition is clearly observed at 7 = j c = 1/2. If we were to use 
the critical point (j) c = 7r, the curve describing the quadrature of the radiation field would 
change sign. 
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FIG. 1: (Color Online.) Critical values (lower, blue) and 9 C (upper, red) as functions of the 
interaction strength 7, for uja — 1- 

From the relations in Eq.Q we obtain, in the superradiant region, 



shifl c , 

; COS (b c . ( 5 ) 

'N v "a/2 cos# c y v ; 

This is plotted in Fig. [2} Notice that the normal regime is described only by the origin 
(q c = 0, 6 C — 0), and that the curve is valid for all \ j\ > j c (as this parameter drops out 
from Eq.([5])). It is also valid for any number of atoms. As we will show below, the same curve 
is obtained for symmetry-adapted variational states as well as for exact quantum solutions; 
in this sense the curve is universal. Furthermore, the first-excited states fall along the same 
path. 

The Hamiltonian ([!]) shows a parity symmetry given by 



[e^ A , = 



(6) 



with A = a^a + J z + ^JJ 2 + \ — |, the total excitation number operator. This allows for the 
classification of its eigenstates in terms of the parity of the eigenvalues A = v + m + j of 
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FIG. 2: Universal curve q c /yN versus # c , for uoa — 1- 

A, where m, and j(j + 1) are eigenvalues of a^a, J Z) and J 2 , respectively. This allows us 
to build, through the projection operators P± — \ (i±e^ 7rA ), symmetry- adapted coherent 
states (SAS) that preserve the symmetry of the Dicke Hamiltonian, i.e., have a definite 
parity of the total number of excitations. These are given by [7] 

\a, C)± = A4 ® |C> ± | — «> (8) | - 
where the normalization factors A/± are 

i-ICI 2 







(7) 



•^" 2 = 2 ^l±exp(-2H 2 ) Vi + K|2 

The energy surface associated to SAS is given by [7j 
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7V N 



(8) 



= 4(p 2 + . 2 ){i- t 



<{ (cos 6*) 



±i 



±e ± (f 2 +9 2 )(cos^)^ 
tan 2 6> cos 6 



+ ^27V 
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±p tan 6* sin (f> + q e p2+q2 sin 9 cos (cos #) -7V 



(9) 



e^ 2 (cos^)-^±l 

It is straightforward to show that in the thermodynamic limit both surfaces {H)+ and (H)_ 
reduce to Eq.Q. 

For a given number of particles N and a fixed coupling parameter 7, we determine the 
values of q C) p C) 6 C) (j) C) where the energy surface for the SAS states has its minimum value, 
in the same manner as done for the coherent states (CS). The minima of the SAS energy 
surface occurs for p c = and (j) c = 0, n. The values of q c and # c , which now depend strongly 
on N : are found numerically and they are plotted as a function of 7 in Fig. [3] for the even 
and odd variational states. 
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FIG. 3: (Color online.) Critical values of q/y/N (lower, blue), and 9 (upper, red) as functions of 
the interaction strength 7, for <j) c — 0, uja — 1 and N = 20 atoms. Left: even SAS (ground state); 
right: odd SAS (first excited state). 

The discontinuity in the phase space variable for the even states marks the phase tran- 
sition. This phase transition arises from a jump of the global minimum in the even energy 
surface (H)+ between two local minima, as 7 varies; this can be seen in Fig. [1] where, for 
three values of 7 crossing the discontinuity, the contour energy levels are plotted as functions 
of q and 0, for (j) c = 0. The jump occurs at 7 = 0.553 for TV = 20. For the odd state, the 
appearance of a phase transition manifests itself through a change in the concavity of the 
associated curves for q c and 9 C . 




FIG. 4: Contours for the even energy surface showing the jump of the global minimum as 

the interaction strength 7 crosses its critical value at the phase transition. Plotted are the critical 
values of q and 0, for N = 20 particles and uja — 1. {E m i n , q m i n ^ Omin} are, from left to right, 
{-10.1887,-0.935972,0.29446}, {-10.1963,-0.964931,0.30329} and {-10.2559,-2.0,0.615064}. 

It is useful to compare the minimum points q c associated to the energy surfaces of the 
CS and the even SAS (cf. Fig. [5]). The quadrature q c shows a discontinuity for the latter 
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at the phase transition, as a function of the coupling parameter 7. This translates into a 
forbidden range of values for 9 C (0.3 < 9 < 0.6) when q c is plotted against it. Note that, 
in the allowed ^-region, the functions q c (9 c ) associated to the energy surfaces of the CS and 
even SAS fall along the same curve. The odd SAS state is also plotted, falls along the same 
curve, has no forbidden regions, and is indistinguishable from that of the coherent state. 
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FIG. 5: Color online. Left: q c /y/~N as a function of 7 for the coherent (continuous, cyan) and 
even SAS (discontinuous, blue) estimations of the ground state. Beyond the phase transition both 
curves coincide. Right: q c /y/N as a function of 9 C for the CS and odd SAS states (continuous, 
green) and even SAS state (discontinuous, blue). 

Since the q c vs. 9 C seems to show a universal behavior, it is convenient to study the exact 
quantum solution. We propose the following correspondence between the critical (g c , 9 C ) 
parameters of the classical phase space and the quantum operators: 

q c — ► TVP>) #c — ► arccos ( Z y^) , (10) 

where j = N/2 and the =p corresponds to <j) c = 0, n. As the system is not integrable, 
we solve numerically the Hamiltonian eigenvalue equation for the ground- and first-excited 
states, separating the even- and odd-parity cases. 

Fig.[6]shows the results. The figure on the left shows, in the vicinity of 7 C , the dependence 
with 7 of the quadrature (or expectation value of the number operator), for A = 20 (lower 
continuous, gray dotted line) and N = 60 (upper continuous, black dotted line). These 
are compared with those for the SAS state at the same values of N (discontinuous, dotted 
blue lines) and with the result for the coherent state (upper continuous cyan line). Note 
that the quantum and SAS solutions agree perfectly at (N p h) = q c = 0, and approximate 
very well beyond the phase transition. While the SAS plot is discontinuous at the phase 
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FIG. 6: Left: Dependence with 7 of the quadrature (or expectation value of the number operator), 
for N = 20 (lower continuous, gray dotted line) and N = 60 (upper continuous, black dotted line), 
compared with those for the SAS state at the same values of N (discontinuous, dotted blue lines) 
and with the result for the coherent state (upper continuous cyan line). Right: Quantum state q c 
vs. 9 C (a/ (aJa)/j vs. arccos(( J z )/j)) behavior, falling exactly on the universal curve together with 
the CS and SAS states. 

transition, the quantum solution shows an inflection point. The CS agrees beyond 7 = 7 C , 
but fails to reproduce the results well enough in the normal region. As TV increases, j c —> 0.5 
(the value for the coherent state), and the jump in the SAS curves becomes smaller, the 
curves approaching more and more its quantum counterparts. The figure on the right shows 
the quantum state q c vs. 9 C (yj (aJa)/j vs. arccos((— J z )/j)) behavior for N = 20, falling 
indistiguishably on the universal curve together with the CS and SAS states. 

The quantum phase transition was found through the fidelity between the state at 7 + £7 
and the state at 7 [12] 

r=\Mi) l#y + *r)> I 2 , (ii) 

which at the transition acquires its minimum value, by taking £7 = 0.001; for N = 20 we 
find the transition at 7 C = 0.567 ± 0.001, comparing very well with that obtained for the 
SAS states, even for this small N (cf. Fig. [7]). Another method is to use its second derivative, 
the so-called fidelity susceptibility, which, besides pinpointing the phase transition, it allows 
the determination of the scaling behavior with respect to the number of particles [T2HT1] . 
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FIG. 7: Quantum phase transition determined through the fidelity of 7- neighbouring states. For 
N = 20 and uja — 1 the quantum phase transition occurs at 7 = 0.567 ± 0.001. 

III. CONCLUSIONS 

A universal curve for all values of the coupling parameter 7 and any number of atoms TV 
is obtained for the phase space variables q versus 0, valid for the ground- and first-excited 
states. The universal character is irrespective of whether one uses coherent CS states, 
symmetry-adapted SAS states, or the exact quantum solution, even though these quantities 
are arrived at via very different methods. The gap in this curve present for the SAS states 
closes in the thermodynamic limit. 

The critical value 7 C at which the phase transition occurs is calculated, for the exact 
solution, by using the fidelity between neighboring states in the coupling strength. This 
agrees very well, though not identically, with the critical value given by the SAS states. 
For instance, in the case = 20 we find j c \sas = 0.553 while jdquant = 0.567. These 
values become even closer as N increases. This will be studied thoroughly by means of the 
fidelity susceptibility elsewhere, together with the properties of the ground- and first-excited 
states at the phase transition, for a finite number of atoms. These results are particularly 
relevant in atomic physics and quantum optics using superconducting g-bits [15j, where 
having semi-analytical solutions for small values of N is very useful. 

While for the coherent states, or mean-field procedure, the phase transition is determined 
properly only in the thermodynamic limit, for our SAS variational states the behavior of 
both, the ground- and first-excited states, are determined by means of jumps in the phase 
space variables as functions of the coupling parameter between the electromagnetic field 
and the matter, and may be calculated for a finite number of particles. As expected, in the 
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thermodynamic limit all solutions converge to the mean field result. 



Acknowledgements 

This work was supported by CONACyT (under project 101541) and DGAPA-UNAM 
(under project IN102811). 



[1] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[2] K. Hepp and E. H. Lieb, Phys. Rev. A 8, 2517 (1973). 

[3] K. Baumann, C. Guerlin, F. Brennecke, and T. Esslinger, Nature 464, 1301 (2010). 
[4] D. Nagy, G. Konya, G. Szirmai, and P. Domokos, Phys. Rev. Lett. 104, 130401 (2010) 
[5] B.M. Garraway, Phil. Trans. R. Soc. A 369, 1137 (2011). 

[6] O. Castanos, E. Nahmad-Achar, R. Lopez-Peria, and J.G. Hirsch, Phys. Rev. A Rapid Comm. 
83, 051601(R) (2011). 

[7] O. Castanos, E. Nahmad-Achar, R. Lopez-Peha, and J.G. Hirsch, Phys. Rev. A 84, 013819 
(2011). 

[8] K.T. Hecht, The Vector Coherent State Method and Its Applications to Problems of Higher 

Symmetries (Springer, Berlin 1987). 
[9] FT. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, Phys. Rev. A6, 2211 (1972). 
[10] O. Castanos, E. Nahmad-Achar, R. Lopez-Peha, and J.G. Hirsch in Symmetries in Nature 
AIP Conference Proceedings 1323, p. 40-59 Eds. L. Benet, P.O. Hess, J.M. Torres, K.B. Wolf 
(Melville, New York, 2010). 
[11] R. Gilmore, Catastrophe Theory for scientists and engineers, (Wiley, New York, 1981). 
[12] P. Zanardi and N. Paunkovic, Phys. Rev. E 74, 031123 (2006). 
[13] S-J. Gu, Int. J. Mod. Phys. B 24, 4371 (2010). 

[14] O. Castanos, R. Lopez-Peha, E. Nahmad-Achar, and J.G. Hirsch in XXXV Symposium on 
Nuclear Physics IOP Journal of Physics: Conference Series (2012, to be published). 

[15] J.Q. You and F. Nori, Nature 474, 589 (2011); J. Clarke and F.K. Wilhelm, Nature 453, 1031 
(2008). 



10 



